The final mass and spin of black hole mergers 
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We consider black holes resulting from binary black hole mergers. By fitting to numerical results 
we construct analytic formulas that predict the mass and spin of the final black hole. Our formulas 
are valid for arbitrary initial spins and mass ratios and agree well with available numerical simu- 
lations. We use our spin formula in the context of two common merger scenarios for supermassive 
galactic black holes. We consider the case of isotropically distributed initial spin orientations (when 
no surrounding matter is present) and also the case when matter closely aligns the spins with the 
orbital angular momentum. The spin magnitude of black holes resulting from successive generations 
of mergers (with symmetric mass ratio n) has a mean of 1.7377 + 0.28 in the isotropic case and 0.94 
for the closely aligned case. 

PACS numbers: 04.25.D-, 04.25.dg, 04.70.-s, 97.60.Lf, 98.62.Js, 98.65.Fz 



Introduction. Barring the influence of surrounding 
matter or third objects, two gravitationally bound black 
holes (BHs) (with masses M a , Mb and spins S a ,Sb) will 
orbit around their common center of mass emitting gravi- 
tational radiation, which carries away energy, momentum 
and angular momentum. This radiation will circularize 
the orbits of the progenitor BHs and eventually shrink 
the orbits until the BHs merge and form a BH of mass 
Mf and spin Sf. Thus, the initial state is described by 
eight parameters: the mass ratio q = Mb/M a , the dimen- 
sionless spins a — Sa/M 2 , b — Sb/M 2 of the initial BHs 
and the dimensionless angular velocity u> = 0(M a + Mb). 
Here w specifies which point on the possibly very long in- 
spiral trajectory is used as initial point. After the merger 
the final BH is characterized by seven parameters, the fi- 
nal mass m = M f /(M a + M b ), the spin s = S f /M 2 f and 

the kick velocity k. Predicting the final m and s from 
the initial parameters is of great importance in many 
astrophysical merger scenarios. Boyle, Kesden and Nis- 
sanke |l|, l2j propose to describe any of the final param- 
eters as a Taylor expansion in the six initial spin com- 
ponents. These spin expansions (with coefficients fitted 
to date from numerical simulations) are based on the as- 
sumption that any dimensionless final quantity must be 
a function of the eight initial parameters. In this paper 
we use their expansions for the final dimensionless mass 
and spin. Note, however, that all expansion coefficients 
still depend on q and w, since we Taylor expand only in 
cti and hi. I.e. even if we know the coefficients for a par- 
ticular mass ratio and a particular initial angular velocity 
we cannot predict any final quantities for different mass 
ratios or initial velocities. At first glance this seems to 
severely limit the usefulness of these expansions. How- 
ever, as we will see, we have come up with a particular fit 
for the g-dependence, which seems to work rather well. 
Also, we will show that the final spin magnitude depends 
only weakly on 10. The orientation of the final spin, how- 
ever, does depend on oj as one would expect due to spin 
precession of the individual BH spins if one starts at a 
different to. Nevertheless, our approach can even give ap- 



proximate spin orientations if we start with values for uj 
like the ones used here, which are typical for the current 
state-of-the-art simulations performed by most groups. 

We use our formulas to calculate the probability den- 
sity of the magnitude of the final spin for successive gen- 
erations of BH mergers. We consider "gas-rich" or "wet" 
mergers, where a circumbinary disk surrounds the binary, 
and also "gas-poor" or "dry" mergers, where no matter 
is present [3, [4] . 

Spin expansions. Our coordinate system is such that 
the z axis is perpendicular to the initial orbital plane. 
The center of mass is initially at rest at the origin, with 
the x axis along direction of the momentum of BH b. 
The y axis is along the line connecting the BHs with BH 
a located at y > 0. 

To construct formulas aimed at predicting the final BH 
mass and spin, we follow the method described in [2j . To 
linear order, the spin expansions for the final mass and 
spin are: 
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m u + {m ai a z + m bi b z ) 



{sfa x + 8?b x ) + (s^ay + s b x %) 
{sfa x + sfbz) + (sfa v + sfby) 
sl + {sfa z + s b z l b z ), 



(1) 



where the coefficients are functions of the mass ratio q 
and the initial orbital angular velocity u>. Notice that 
by using symmetries such as parity or exchange many 
terms that would appear in an unconstrained Taylor ex- 
pansion have been dropped. In addition, all the co- 
efficients enclosed in common brackets are related, by 
m bl (q,u>) = m al (l/q,ui) and s^ (q,u) — s^ 3 (l/q,uj). 

We ignore the dependence on u> for the time being and 
focus first on the equal mass case. Note that for q = 1 
the above mentioned relation between coefficients implies 
that all terms enclosed in common brackets have equal 
coefficients. In order to determine these coefficients we 
have performed 10 numerical simulations of equal mass 
binaries with spins of magnitudes between 0.1 and 0.27 
with orientations as in Table XI of [2|. All 10 simulations 



-0.013 ±0.007 
sf = sf = 0.187 ± 0.002 = +sf = +sf 
sf = sf = 0.028 ± 0.002 = -sf = -sf 
s° z = 0.686 ± 0.004 sf = s" =0.15 ±0.03 

TABLE I: Equal mass coefficients up to linear order. 



start with u) = 0.05 and are performed using the "mov- 
ing punctures" method @, @| with the BAM code 0, Q 
which allows us to use moving nested refinement boxes. 
We use 10 levels of 2:1 refinements. The outer bound- 
aries are located 436M away from the initial center of 
mass (M being the sum of the initial BH masses), and 
our resolution ranges between 8M on the outermost box 
to M/64 near the BHs. Apart from the use of sixth or- 
der stencils in the interior of the boxes [9|, our setup 
and methods to determine spin and mass are very sim- 
ilar to the simulations reported in [10L 111 , 



12| . Our 10 



simulations are in principle sufficient to determine all co- 
efficients appearing in Sj and m up to quadratic order. 
However, since the numerical errors in determining the 
final mass and spin are 0.1% and 0.5% respectively, and 
since all quadratic coefficients are small, the errors in all 
quadratic coefficients are more than 100%. For this rea- 
son only the linear coefficients are listed in Table [TJ In 
order to verify that linear expansions with the coefficients 
in Table UJ give reasonable results we have performed 10 
more equal mass runs with spins of magnitude 0.75 with 
arbitrary orientations. As one can see from the first 10 
lines in Table [Til the agreement between our additional 
runs and the values predicted by the expansions is quite 
good. This demonstrates that for q = 1 we can trust the 
expansions with our coefficients. As we can see from Ta- 
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s^ , which is not predicted 



by the symmetries used to derive the spin expansions. 
Thus, in terms of the initial spins S a , Sb and final spin 
Sf the x and y components of Eq. {TJ) can be rewritten 



S f,x 

S f,y 



aR 






Sb,y 



(2) 



with 



R = 



cos/3 
sin/3 



— sin/3 
cos/3 



/3 = 0.15, a = 0.685. (3) 



This means that the final spin components in the xy- 
plane are given by rotating and scaling the initial spins. 
The scale factor is a — 0.685 and the rotation matrix 
is R. This result is physically reasonable, because we 
know that the initial spins will precess and also radiate 
some angular momentum. If Eq. ^ also held for unequal 
masses it would imply a certain form for the mass ratio 
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TABLE II: Test simulations: the columns show the initial 
mass ratio q, the progenitor spins a* and bi, the final spin 
magnitude s and mass m, and our predictions s v and m v . 
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However, the latter cannot be true for extreme mass ra- 
tios, since than we would not get s = a for q = 0. In order 
to get the extreme mass ratio limit correctly, we modify 
the mass ratio dependence somewhat and assume that 
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[g(q,c)a x +g(l/q,c)b x 
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(c+l) 2 
{c+qf 



1 



= 0.762, 

(6) 

are chosen such that s '= a for q = and s = b for q = oo. 

So far we have only considered the components in xy- 

plane. For equal masses sf = sf so that s z in Eq. (JTJ) 

simplifies. As above we introduce a (/-dependence given 



by 



= s Q z {Awr 1 +16{l-w)r ] 2 )+ s f [g(q,c 3 )a z + g(l/q,c 3 )b z } , 

(7) 
= 0.632. The r\ dependence 



where c 3 = 

of the leading coefficient is inspired by post-Newtonian 
(PN) expressions, and has been shown to fit runs without 
initial spins very well [13(. The constant w = 1.26 is fitted 
using the results for unequal masses discussed below. 

Similarly, we build a mass ratio dependence into the 
final mass formula and write 

m = 1 + (m° - 1)4:7] + m al 16r] 2 (a z + b z ) . (8) 

The leading term is inspired by the fact that the binding 
energy is proportional to 77, and the linear term is chosen 
such that it has the correct limit in the extreme mass 
ratio case. 

While the q dependence we have introduced so far re- 
produces our fits for the equal mass case and simultane- 
ously gives the correct answers in the extreme mass ratio 
cases, it is not clear how well our formulas perform for 
intermediate mass ratios. For this reason we have per- 
formed 21 more simulations with q different from unity. 
As one can see from Table ITT1 the predictions are at most 
0.4% from the numerical results. 

So far we have only kept linear terms in our expansions. 
We find that our formulas agree with a wide range of 
test runs. Yet, for initial spin magnitudes close to 1 our 
formulas deviate from the extrapolated values l_2j, |l4j, |l5j 
for the minimum and maximum possible final spin. These 
problems can be fixed if we add some quadratic terms to 
s, alone: 



- Jo) 



(4wrj + 16(1 ~w)r) 2 ) 



+sf [g{q,c 3 )a z + g(l/q,c 3 )b z ] 

+16kr 1 2 [(a x + b x f + (a y + b y f - (a z + b z ) 2 }(9) 

These quadratic terms approximately reproduce most of 
the very small and very uncertain quadratic coefficients 
in s z . The coefficient k = 0.008 is chosen such that we 
get the best overall agreement with our numerical simu- 
lations. 

Results. Our particular mapping of initial masses and 
spins into the final mass and spin given by Eqs. HJ, |5]), 
([5]) and ^ was fitted for the initial orbital angular veloc- 
ity uj = 0.05. If we start with the same initial spin com- 
ponents but at lower initial angular velocity (i.e. larger 
separation) the individual spins and the orbital plane will 
have precessed by the time we reach uj = 0.05. In this 
case we cannot expect that our formulas will predict the 
final spin components if we simply use the initial spin 
components. However, the final spin magnitude should 
still be approximately correct, since PN calculations [16j 
demonstrate that the spin magnitudes are conserved at 
2PN order. This expectation is borne out for the follow- 
ing additional q = 1 test run. It starts with uj = 0.03, 



a = (-0.637, -0.226, 0.325), b = (-0.517, -0.543, 0.025), 
and yields a final spin s = (-0.226,-0.146,0.746) with 
magnitude s — 0.793. Our fitting formula predicts 
?= (-0.194,-0.176,0.753) with magnitude s = 0.797. 
This shows that the predicted magnitude is correct up to 
an error of 0.6%, while the components only agree if the 
predicted vector is rotated. This rotation comes from the 
precession during the time it takes to go from uj = 0.03 
to uj = 0.05, during which the system completes about 4 
orbits. 

We have also compared with the numerical results pub- 
lished in [TIL E3, EE [IE E3, EE [S H3] , which all start 
from an uj that is not too far from 0.05. The average de- 
viation of all these results from our predictions is about 
1%. Hence our formulas gives useful predictions for the 
final mass and even for the spin components if one starts 
near uj = 0.05. However, if one is ever interested in the 
final spin components for a much larger initial separation 
with uj <C 0.05, one can still evolve the spins using PN 
theory up to uj = 0.05 and then use our formula. This 
eliminates the need for expensive numerical simulations. 

Other groups have presented formulas that, like ours, 
attempt to predict the final spin of the merger. The an- 



alytic estimate of Buonanno et al. [2l| can give the final 
spin magnitude to within few percent with larger devi- 
ations for spins close to anti alignment. Starting from 
a number of assumptions, Rezzolla et al. [17| developed 
a more accurate formula with coefficients fitted to nu- 
merical results. One of their assumptions is that the 
components of the final spin in the initial orbital plane 
are obtained by summing the components of the initial 
spins i.e., Sf iX / y — S a _ x / y + St, yX / y . In our simulations 
we observe that this assumption is violated. From the 
discussion around Eqs. [2] and [3] it is clear that these com- 
ponents get slightly rotated and shortened by a factor 
a = 0.685 during the merger. For instance, the run from 
the first line in Table HT1 we find a value of s xy = 0.291 for 
the in-plane component of the spin. Our formula ([5]) pre- 
dicts a close value of s XV:P — 0.284, the approach in [Uj 
leads to s xy ^ = 0.375 which is about 30% too large, while 
both approaches predict almost the same spin magnitude. 
Thus our formulas improve the final spin orientation. 

Black hole mergers. Using Eqs. ([5]), (O and © we 
can study the properties of the spin of BHs produced by 
successive binary BH mergers. Two types of scenarios 
are likely to witness mergers: one in which the two BHs 
carry with them dense surrounding matter ( "wet" merg- 
ers) and another where the progenitors meet in relatively 
empty space ("dry" mergers) 3]. In the former case the 
progenitors spins are likely to reach the merger aligned 
with the orbital angular momentum 22j, while in the 



latter the same spins are bound to be isotropically ori- 
ented [3j. In our models we assume that the initial spin 
directions either have a uniform probability distribution 
("dry" mergers) or that the longitudinal angle 9 obeys a 
normal distribution centered in the direction of the or- 




FIG. 1: Probability density (for different mass ratios) of the 
final BH spin magnitude after 4 generations of "dry" mergers 
(top) and "wet" mergers (bottom). 
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FIG. 2: Mean value of the final spin magnitude as a function 
of the symmetric mass ratio r\. The vertical bars are the 
corresponding standard deviations. 



bital angular momentum ("wet" mergers) [23|. For the 
spin magnitude we assume that the original progenitors' 
spin magnitude is uniformly distributed in the interval 
[0,1]. The probability density for the final spin mag- 
nitude after merger can then be obtained with Monte 
Carlo simulations. Note that the u> dependence of our 
formula is removed by the integration over the spin ori- 
entation. This procedure can be iterated for successive 
generations assuming that the magnitude of the progen- 
itors spins has the probability density of the previous 
generation. The resulting probability distribution con- 
verges quickly to the curves shown in Fig. [1] for the mass 
ratios q = 0.4, 0.6, 0.8 and 1.0. The mean values for the 
probability densities as a function of the symmetric mass 
ratio r\ are shown in Fig. [5] with the vertical bars corre- 
sponding to the standard deviations. The straight lines 
are linear fits of the data points. We find that successive 
"dry" mergers produce final spins with a mean value of 
(s)Dry = 1.73?7 + 0.28. For q = 1 the final spin does not 
exceed s = 0.954 (as our formula predicts for maximal 
aligned spins). Successive "wet" mergers produce a spin 
around 0.94, with very little spread. Similar results have 
also been found by Berti et al. [4] using the spin formula 
in [171 ] . This agreement is expected since our spin formu- 
las give magnitudes in close agreement with [171]. The 
main difference is in the spin orientation which has been 
integrated out. 



Discussion. By fitting to numerical results we con- 
struct formulas [Eqs. (J8]), ([5]), dH)] that predict the mass 
and spin of the final BH coming from binary BH mergers. 
We use them to determine the probability distribution of 
the final spin magnitude (Figs. [TJ \2§ after several gener- 
ations of mergers of either "dry" or "wet" mergers. 
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